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1 Introduction 

In recent years there has been tremendous progress in constructing supersymmetric lattice 
gauge theories, as reviewed by ref. [1] and references therein. In particular, a lattice 
formulation of four-dimensional J\T = 4 supersymmetric Yang-Mills (SYM) theory has 
been constructed, which exactly preserves a single twisted supercharge Q at non-zero lattice 
spacing [2-12 ]. 1 Initial studies, both analytical and numerical, have shown that this lattice 
theory possesses several remarkable properties: 

• The renormalization of the theory is highly constrained by gauge invariance, the exact 
Q supersymmetry and a global £5 symmetry. These symmetries ensure that at most a 
single marginal coupling may need to be tuned to correctly recover all 16 supercharges 
and the SU(4) R symmetry of AT = 4 SYM in the continuum limit [ 8 , 10]. 

• The moduli space of the lattice theory is the same as that of the continuum theory, 
and is not lifted to any order of lattice perturbation theory [ 6 ] . 

• The /3 function vanishes at one loop in lattice perturbation theory [ 6 ]. 

1 Other approaches to studying J\[ = 4 SYM numerically include refs. [13-19]. 
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• Although the pfaffian that results from integrating over the fermion fields in the 
lattice path integral is potentially complex, pfaffian-phase-quenched computations 
display no indication of a sign problem [9]. This allows numerical computations to 
employ conventional Monte Carlo algorithms [11]. 

• The static potential of the lattice theory is coulombic at both weak and strong cou¬ 
pling, in agreement with continuum expectations [7, 9, 12]. 

Given the central role J\f = 4 SYM plays in the AdS/CFT correspondence that relates it 
to quantum gravity, it is important to carry out large-scale lattice calculations to investigate 
the theory away from the regime of weak coupling and for arbitrary numbers of colors N. 
Despite the encouraging results described above, the formal action constructed in refs. [2, 4] 
possesses certain features that obstruct such numerical calculations. The primary issue, 
explored in ref. [9], is that this lattice system possesses flat directions that destabilize the 
vacuum of the theory. In particular, the theory necessarily involves U(AQ = SU(A^) <8> U(l) 
gauge invariance, and flat directions associated with the U(l) sector produce especially 
severe lattice artifacts. This behavior may be understood by noting that constant U(l) 
shifts of the fields A'(x) —>• X(x) + cl leave the action S unchanged for any x-independent 
field configurations X{x) even when S'fY] ^ 0. This invariance should be contrasted with 
that of the usual SU(iV) flat directions that correspond to constant shifts by elements of 
the Cartan subalgebra on x-independent field configurations for which the action S = 0. 
In other words, the SU(iV’) flat directions involve only supersymmetric vacuum solutions 
while the U(l) flat directions can play a more general role. Of course the continuum theory 
admits a trivial truncation to SU(-ZV), but this is not true for the lattice theory at non-zero 
lattice spacing. 

In ref. [9] the fluctuations of these U(l) modes were controlled in a way that violated 
the exact supersymmetry. Although the resulting soft supersymmetry breaking was under 
control for ’t Hooft couplings Aj a t < 5, and suppressed oc 1 /N 2 , it limited the exploration 
of stronger couplings and produced a slow approach to the continuum limit [12]. In this 
paper we show how these problems may be solved through a construction that regulates 
the U(l) modes in a manner compatible with the Q supersymmetry. The basic idea is to 
modify the equations of motion of an auxiliary field, which determine the moduli space of 
the system. By maintaining the Q supersymmetry we ensure that Q-invariant observables 
such as the partition function can still be computed by restricting the functional integrals 
to a moduli space that has been modified to include only the SU(A r ) sector. 

We begin in the next section with a brief review of the difficulties faced by the formal 
supersymmetric lattice action, summarizing how these were addressed by the unimproved 
action used in the past, at the cost of soft supersymmetry breaking. In section 3 we describe 
the general framework of Q-invariant modifications to the moduli equations, which we use 
to construct our improved lattice action. We test this improved action in section 4, finding 
that it dramatically reduces violations of Ward identities on small L 4 lattice volumes and 
exhibits much more rapid approach to the continuum limit as L increases. Although the 
improved action is almost twice as computationally expensive as the unimproved action, 
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these benefits are so substantial that we recommend using the improved action in future 
studies of lattice Af = 4 SYM. 

One may also consider alternate modifications of the moduli equations, for which this 
procedure produces other possible lattice actions. We discuss one interesting example in 
appendix A, which turns out not to behave as well as the improved action we focus on. 
Although further explorations in this direction may be worthwhile, the improved action 
we present here appears hard to beat in the context of practical numerical calculations. 
In particular, in appendix B we argue that the symmetries of the lattice system forbid 
all dimension-5 operators, which allows the new action to be effectively 0(a) improved. 
Finally, in appendix C we exploit the improved action to provide a brief update on the 
absence of a sign problem in lattice J\T = 4 SYM. 

In the course of this paper we will discuss several different lattice actions, which we 
briefly define here to avoid confusion: 

• The formal action <Sf orma i directly discretizes continuum twisted AT = 4 SYM, exactly 
preserving the Q supersymmetry but suffering from difficulties summarized in the 
next section. 

• The unimproved action S un i mp used in refs. [9, 11, 12] stabilizes numerical calcula¬ 
tions through two new terms, each of which softly breaks supersymmetry. 

• The improved action <Si mp introduced in this work replaces one of the two soft Q- 
breaking terms of the unimproved action by a supersymmetric construction, removing 
the dominant source of supersymmetry breaking. 

• The over-constrained action S over discussed in appendix A implements all necessary 
stabilizations in a Q-invariant manner, but this turns out to introduce new difficulties. 


2 Difficulties of the formal supersymmetric lattice action 

The starting point for lattice Af = 4 SYM is the direct and exactly Q-supersymmetric 
discretization of the continuum twisted action, 


•Sformal Sexact T S'closed 


5 , 


exact 


N 

2Alat 

N 




Q (Xabin)V^U b (n) + ri(n)v[ W a (n) - \i](n)d{n) 


^closed — ^ x E Tr I ^abcde Xde (bl T /Tq T f^b T /^c)T* c Xab(j^) 
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( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


with repeated indices summed. All fields X = Y2a =i T A X A are expanded on a basis of 
anti-hermitian generators of the group, Tr [t a T b ] = —S AB . Here U a with a = 0 • • • 4 are 
complexified gauge links that also contain the scalar fields and d is a bosonic auxiliary field, 
while r/ and Xab = — Xba account for 11 of the 16 twisted fermion field components. The 
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lattice finite-difference operators are [4, 5] 

V^f b (n) = U a (n)f b (n + J2 a ) - f b (n)U a (n + J2 b ) 

} /a(n) = f a (n)U a (n) - U a (n - jj, a )fa(n - fi a ) (2.4) 

£*1 +) Mn) =U a {n + fib)fb( n ) - f b (n + /2, a )U a {n ) 

T>i ) fab(n) = f ab (n + ]u, c )U c {n ) - U c {n + j2 a + ]u b )f ab {n). 

This lattice action is invariant under the Q supersymmetry 

Q Ua(n) = i/ja(n) Q V> a(n ) = 0 

Q Xab{n) = -Fabin) Q U a {n) = 0 (2.5) 

Q r)(n ) = d(n) Q d(n ) = 0 


where i/j a are the remaining 5 fermion field components and the complexified field strength 
is T ab = V^U b = — F ba ■ From eq. 2.5 we see Q 2 = 0, which establishes the invariance 
of iSexact - The invariance of S^iosed follows from an exact lattice Bianchi identity. If we 
integrate out the bosonic auxiliary field d through its equation of motion d = T>" a ' > 1A 0 , 
eq. 2.2 becomes 


'exact — 


N 

2Alat 


E Ti - 


T‘ a fe(n)J ab (n) - Xab{n)V { ^ ] ^(n) - rj(n)V { a VaW ( 2 -6) 

+ \(vVu a {n )) 2 . 


This Q-invariant action possesses certain features that interfere with numerical calcu¬ 
lations. These features are direct consequences of the structure of the continuum theory 
and the requirement of exact Q supersymmetry at non-zero lattice spacing. First note that 
eq. 2.5 forces the complexified gauge link variables U a to reside in the algebra gl(N,C) of 
the group U (N), in correspondence with i/j a . 2 Therefore targeting the correct continuum 
theory requires that the links have the expansion U b (n ) = + A b {x) + 0{a) in some 

suitable gauge, where a is the lattice spacing and A b is the complexified gauge field of 
the twisted continuum theory. We can interpret the unit matrix in this expansion as the 
vacuum expectation value (vev) of an imaginary U(l) mode, implying that we must lift 
the corresponding flat directions. To do so we need to add to the action a scalar potential, 
typically (^Tr \U a Ua\ — l) , which will also regulate the usual SU(1V) flat directions. 

However, such a scalar potential does not affect real U(l) modes, which cancel out 
of U a U a . While these U(l) modes simply decouple in the continuum theory, in ref. [9] we 
showed that on the lattice they lead to unphysical confinement at strong coupling, indicated 
by the condensation of lattice-artifact U(l) monopoles. We therefore need to add another 
regulator to the action, which we take to be some function of {detV ab ~ 1), where 


V ab (n) = U b {n)U a (n + fi b )U b (n + 'p, a )U a {n) (2.7) 

2 The lattice path integral over gauge fields that live in the algebra involves the flat measure rather than 
the usual Haar measure familiar from lattice QCD. This formulation remains gauge invariant since the 
Jacobian resulting from a gauge transformation on DU cancels against the corresponding quantity for DU. 
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is the oriented plaquette in the a-b plane. Since this term does not affect the SU(JV) flat 
directions, we still need to retain the scalar potential to fully stabilize numerical calcula¬ 
tions. 

Finally, the fermion action possesses an exact U(l) zero mode corresponding to the 
shift symmetry 77 —> 77 + cl at, with c a constant Grassmann parameter. We typically lift 
this zero mode by imposing anti-periodic (thermal) temporal boundary conditions (BCs) 
for the fermion fields, which explicitly breaks the Q supersymmetry since the bosonic fields 
are always periodic in all directions. This source of supersymmetry breaking is controlled 
by the temporal extent of the lattice, Nt- Previous studies observed negligible effects even 
for Nt = 4 [9], though this may be a non-universal feature of the unimproved action then 
employed. We will revisit this issue for the improved action in section 4. 

In summary the formal lattice theory suffers from problems that have their origin in 
a U(l) sector that necessarily arises in the lattice construction and that does not decouple 
from the physical SU(IV) sector at non-zero lattice spacing. This U(l) sector possesses 
flat directions that can destabilize lattice calculations. To avoid instabilities and ensure 
that the lattice theory possesses the correct continuum limit we must regulate these U(l) 
flat directions. In the past, we did so by adding the simplest possible scalar potential and 
plaquette determinant terms directly to eq. 2 . 1 , as 

‘S'unimp = 'S'formal H~ Ssoft 

SWt = (77 Tr pa{.n)u a (n)\ - 1 

Since these terms involve only bosonic and not fermionic fields, they both introduce soft 
supersymmetry breaking that vanishes in a controlled way in the limit ( 7 ,t, k) —> (0, 0) [9, 12]. 
We now describe a new method that allows us to introduce such deformations in a manner 
compatible with the exact Q supersymmetry. 


|detT’ a b(n) - 1 | 


n a<b 


3 Supersymmetric deformations of the moduli space 
3.1 General framework 


The key observation reported by this paper is that one can modify the structure of the 
moduli equations to enforce new vacuum conditions without breaking the Q supersym¬ 
metry. This procedure can be used to regulate the problematic U(l) flat directions. We 
proceed by deforming the Q rjV a 


] Ua term in 5 exa ct (eq. 2 . 2 ), 


Q Tr 


ri(n)(V a U a (n 


Q Tr 


77(77.) 


V { a ] U a (n) + GO(n 


(3.1) 


where G is a new tunable coupling and 0(n) is a gauge-invariant bosonic operator. Al¬ 
though we only recover lattice AT = 4 SYM in the limit G —> 0, the deformed S' exSiCt part 
of the lattice action remains Q exact. 

Performing the Q variation in eq. 3.1 produces 


Q Tr 


77 ( 77 ) } Z4(n) + GO(n 


= Tr 
- Tr 


d{n) [V { a ^lA a (ri) -\- G(D(n 


77(77) ( v [ VaCrc) + GQ 0(n 


(3.2) 
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In combination with the Q r]d = d? term in eq. 2.2 we obtain the modified equations of 
motion 

d{n) = V^U a {n) + GO(n) 1 N . (3.3) 

Since Q rj = d, this deformation also modifies a Q Ward identity, 




Tr [V a U a {n)\ + NG E°(") 


= 1 VO 


E°(’ 


0, (3.4) 


where the first term vanishes due to the sum over the lattice volume and the structure of 
the finite-difference operator, eq. 2.4. After integrating out the auxiliary field, only the 
second line of eq. 2.6 for S^act changes, 


Q> 


exact 


N 

2Alat 


5 > 


Fabir^Fabin) - Xafe(n)P [ ( a f V b ](«) - rj(n)V[ ] ^ a {n) (3.5) 

+ i (V ( a - ] U a (n) + GO(n)I N y - Gr,{n)Q Q{n) . 


3.2 Application to construct the improved action 

Since the plaquette determinant term in 5 so f t (eq. 2.8) produces much larger soft super- 
symmetry breaking than does the scalar potential term [9], we consider 

0(n ) = ^2 (det V a b(n) ~ 1) = 2Re ^ ( detV a b(n) - 1), (3.6) 

a^b a<b 


where the second equality follows from Vi, a = V* b . With this 0(n), the Ward identity in 
eq. 3.4 produces the single constraint 

EE (detV a b(n) — 1) = 0. (3.7) 

n a^b 


Invariance under discrete lattice translations and rotations then implies (Re det V a b(n)) = 
1. Although the deformation of the lattice moduli space does not constrain (Im det V a b(n)), 
the same scalar potential in eq. 2.8 that regulates the SU(1V) flat directions will also address 
this issue, by requiring ^Tr \U a {n)U a {n)\ ~ 1. 

To write down the improved lattice action we need the Q variation of the plaquette 
determinant, which follows from Jacobi’s formula 


d det U a 
da 


[det U n ] Tr 


,_1 dUa 
a da 


(3.8) 


Applying this to the determinant of eq. 2.7, we obtain 


Q [det V a b(n) - 1] = [det V a b(n)\ Tr [U b l {n)^ h (n) +U a 1 (n + ^b)^a{n + /7fe)] (3.9) 
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so that the full improved action becomes 


*Simp — S exact + Closed + 'S’soft 
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(3.10) 
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S’soft = 


N 

2 A lat' 


E E (^f Tr [ u a{n)U a {n)\ - 1 


with S’dosed still given by eq. 2.3. Written in this form, the Q invariance of S' xact is reflected 
by the presence of det V in both the bosonic action and the fermionic action. The plaquette 
determinant terms in the fermionic action, written as Sdet above, break the r/ —> r/ + cIat 
shift symmetry discussed in the previous section, and lift the U(l) fermion zero mode. Sdet 
involves the rj^ a fermion bilinear and can be interpreted as providing an effective fermion 
mass term, with the Tr [//] factor picking out the U(l) sector. 

When implemented in the publicly available parallel software described in ref. [11], Sdet 
also adds 40 inter-node data transfers to the fermion operator, bringing the total to more 
than 100 (many more than lattice QCD with Wilson or staggered fermions). Although 
we no longer need to use anti-periodic temporal BCs to lift the U(l) fermion zero mode, 
we find that these BCs prevent the appearance of a sign problem that results if we use 
fully periodic BCs. Mysteriously, we also measure most observables to be insensitive to 
the BCs and to the presence or absence of a sign problem. We discuss these issues further 
in appendix C and hope to gain further insight from future investigations. As a practical 
matter, the anti-periodic temporal BCs significantly reduce the condition number of the 
fermion operator, accelerating numerical computations. Even so, the additional terms and 
inter-node data transfers cause the improved action to be almost twice as computationally 
expensive as the unimproved action. 

The linear dependence of the new supersymmetric deformation on (det' P a b( n ) ~ 1) 
is another significant difference compared to the [det V a b(n) — 1| 2 used in the unimproved 
action, eq. 2.8. While the latter term was designed to approximately project every plaquette 
to sl(iV, C) at each site of the lattice, the new deformation focuses only on regulating the 
U(l) flat directions. It was not obvious a priori that the linear form would suffice, but 
this does turn out to be the case. However, the linear (det"P — 1) deformation prevents us 
from also including the scalar potential in eq. 3.6 for 0(n). When we do so we find that 
the scalar potential and plaquette determinant largely cancel each other out and fail to 
regulate the flat directions. 

This leads us to retain the same soft supersymmetry breaking scalar potential as 
used in the past. Although this S( oft does much less violence to Q than did the old 
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plaquette determinant term, it is enticing to consider the possibility of obtaining a fully Q- 
invariant stabilized lattice action. In order to do so we must make both the scalar potential 
and plaquette determinant deformations non-negative at each lattice site, for example by 
switching back to the |det V a b — 1| 2 form used in eq. 2.8. We present such a construction 
in appendix A, finding that the resulting action encounters new difficulties and does not 
behave as well as the improved action in eq. 3.10. 


4 Tests of the improved action 
4.1 Lattice phase diagrams 

In this section we present numerical evidence that the improved action does in fact produce 
dramatic improvements compared to the S u nim P used in refs. [9, 11, 12], Each action 
involves two auxiliary couplings: (q, G ) for 5j mp and (/x, k) for S'unimp- Since we only recover 
J\T = 4 SYM upon sending all auxiliary couplings to zero, we wish to carry out calculations 
using the smallest acceptable q, and G / k. If these are made too small, however, the lattice 
calculations will exhibit instabilities, such as the U(l) monopole condensation transition 
discussed in ref. [9] . The minimum acceptable values of the auxiliary couplings depend on 
both the lattice volume L 4 and the’t Hooft coupling Ai at : We must increase q and G / k 
to consider larger Aj a t, but can decrease them as L increases. In addition, for any given 
(L, Ai a t) the minimum acceptable value of G or k depends on /x, and vice versa, leaving us 
with a multi-dimensional parameter space to explore. 

We proceed by roughly mapping out the parameter space using a small 4 4 lattice 
volume, which tends to be more susceptible to instabilities than are lattices with larger 
L. For each ’t Hooft coupling of interest, we scan over small ranges of G for a few values 
of q, based around initial estimates extrapolated from smaller values of Ai a t - For the 
improved action we preferentially reduce /x even if this requires larger G, to minimize soft 
supersymmetry breaking. For the unimproved action we fix k = 0.5 as discussed in ref. [9]. 
Finally, when moving from L = 4 to larger lattice volumes we currently reduce q r \j 1 /L 
while keeping G or n fixed. This scaling for /x can be motivated by interpreting q as a 
mass for the imaginary U(l) modes, and is less aggressive than the /x 2 oc 1/L 4 proposed 
in ref. [9] on the grounds that these modes are constant in space and time. While we have 
found that we can also reduce G on larger volumes, in the next subsection we will see that 
this does not seem to be necessary to rapidly approach the continuum limit where the U(l) 
sector (and thus the plaquette determinant deformation) decouples from the SU (N) target 
theory. 

Further details of the lattice phase diagrams for the improved and unimproved actions 
are not relevant to the present study. Instead, in table 1 we simply list the lattice ensembles 
that produce the results presented in the next subsection, with couplings chosen through 
the procedure summarized above. In this paper we consider mainly gauge group U(2); 
initial explorations for N = 3 and 4 suggest that the same values of /x and G / k should 
be used for any N. 


- 8 - 


Improved action 

Volume Ai a t l l G 

Unimproved action 

Volume Ai a t ^ k 

4 4 

0.5 

0.4 

0.02 

4 4 

0.5 

0.5 

0.5 

4 4 

1.0 

0.4 

0.05 

4 4 

1.0 

0.4 

0.5 

4 4 

2.0 

0.8 

0.05 

4 4 

2.0 

0.5 

0.5 

4 4 

3.0 

0.8 

0.05 

4 4 

3.0 

0.5 

0.5 

4 4 

4.0 

0.8 

0.05 

4 4 

4.0 

0.8 

0.5 

4 4 

5.0 

0.8 

0.10 

4 4 

5.0 

0.8 

0.5 

4 4 

6.0 

0.8 

0.10 

4 4 

6.0 

1.0 

0.5 

4 4 

7.0 

1.0 

0.15 





4 4 

8.0 

1.0 

0.25 





6 4 

1.0 

0.3 

0.05 

6 4 

1.0 

0.3 

0.5 

6 4 

2.0 

0.5 

0.05 





6 4 

3.0 

0.5 

0.05 





6 4 

4.0 

0.5 

0.05 





8 4 

1.0 

0.2 

0.05 

8 4 

1.0 

0.2 

0.5 

8 4 

2.0 

0.4 

0.05 





8 4 

3.0 

0.4 

0.05 





8 4 

4.0 

0.4 

0.05 





12 4 

1.0 

0.15 

0.05 

12 4 

1.0 

0.15 

0.5 

16 4 

1.0 

0.1 

0.05 






Table 1. Auxiliary couplings used to obtain the results in figures 1-7, chosen through the procedure 
summarized in section 4.1. The gauge group is U(2) in all cases. 


4.2 Ward identity violations and continuum extrapolations 

To quantify the severity of supersymmetry breaking, due to both soft breaking terms and 
the anti-periodic temporal BCs for the fermions, we will measure violations of Q Ward 
identities. We begin with the Ward identity resulting from QSformai = 0. Because the 
fermion action is gaussian this Ward identity requires that the bosonic action per lattice 
site take the value (sb) = 9N 2 /2 for gauge group U(1V), independent of the ’t Hooft 
coupling Ai a t. In figure 1 we plot normalized violations of this Ward identity for both 
‘S'unimp and Si m p across a wide range of A[ at on small 4 4 lattices with N = 2. Since we 
have to use larger /r as the coupling increases (table 1) we can see that these violations 
are sensitive to both /r and Ai a t, as expected. The improved action produces roughly an 
order of magnitude reduction in the Ward identity violations, which allows us to reach 
significantly stronger couplings with percent-level Q breaking. 

The anti-periodic fermion BCs produce very little supersymmetry breaking for Simp, 
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Figure 1. Deviations of the bosonic action from its exact supersymmetric value are much larger 
for Sunimp than for Si mp throughout a wide range of’t Hooft coupling Ai a t on 4 4 lattices. Lines 
connect points with fixed p (cf. table 1 ). 



Figure 2. Power-law continuum extrapolations of bosonic action deviations on logarithmic axes for 
5unimp and S'imp with fixed Ai a t = 1- In both cases the deviations vanish oc ( a/L) p in the ( a/L) —> 0 
continuum limit, but the exponent p = 2.42(13) for the improved action is much larger than the 
p = 0.015(4) for the unimproved action. 


which was also the case for the unimproved action [9]. However, the rest of the super- 
symmetry breaking from p ^ 0 is now so small that these BC effects can be relatively 
significant, at least on small 4 4 volumes. For example, switching from anti-periodic to 
periodic BCs for (Ai at ,p,G) = (1,0.4,0.1) roughly halves the normalized bosonic action 
deviations, from 0.0049(8) to 0.0021(5). 

The improvement becomes even more striking when we consider larger lattice volumes, 
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0.02 



Figure 3. Power-law continuum extrapolations of bosonic action deviations on logarithmic axes 
for the improved action with larger 2 < Ai a t < 4 in addition to the same Ai a t = 1 points from 
figure 2. The deviations still vanish oc ( a/L) p at these stronger’t Hooft couplings, with consistent 
p ss 2 listed in table 2. The gap between Ai at = 1 and the other points is due to the larger values 
of n that had to be used for Ai a t > 2 (cf. table 1). 

and in particular if we investigate the approach to the continuum limit as ( a/L) —> 0. In 
figure 2 we plot the deviations of the bosonic action against a/L on logarithmic axes, fixing 
Aiat = 1 for both Sunimp and 5i mp . The right-most results on this plot match the Ai at = 1 
points in figure 1, and as L increases the Ward identity violations decrease for both actions, 
as expected. However the decrease is much more substantial for the improved action, which 
for L = 16 produces violations roughly 500x smaller than those from 5 un i mp . 

Fitting both sets of data in figure 2 to power-law continuum extrapolations oc (a/L) p , 
we find p = 2.42(13) for the improved action. While the Ward identity violations from the 
unimproved action also vanish in the (a/L) —> 0 continuum limit, they approach this limit 
much more slowly, with p = 0.015(4). In figure 3 we check that the improved action behaves 
consistently at stronger couplings, obtaining p « 2 as listed in table 2. The ~(a/L) 2 scaling 
of these results suggests that Si mp is effectively 0(a) improved, i.e. the Q supersymmetry 
is preserved well enough to eliminate all non-negligible effects of dimension-5 operators. In 
appendix B we confirm that exact Q invariance, in combination with the other symmetries 
of lattice J\T = 4 SYM, would forbid all such dimension-5 operators. 

Next, as discussed in section 3, the deformation of the auxiliary field equations of 
motion also introduces a new Q Ward identity for the improved action when G is non¬ 
zero: (1 — Re detR) = 0 from eq. 3.7. Here we average the plaquette determinant over 
all orientations and lattice sites. It is a gauge-invariant quantity associated with the U(l) 
sector, and approaches unity as this U(l) sector decouples from the SU(fV) target theory. 
While the unimproved action does not protect the plaquette determinant through a Ward 
identity, one of the terms in eq. 2.8 suppresses (1 — Re detR) at the cost of soft super- 
symmetry breaking. In figure 4 we plot (1 — Re detR) for S u „\ m r and 5i mp across a wide 
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Figure 4. (1 — Re detR) —>• 0 when the U(l) sector decouples from the SU(iV) target theory. 
Here we average the plaquette determinant over all orientations and lattice sites. This quantity 
is suppressed by one of the soft supersymmetry breaking terms in S'unimp, while S- lmp introduces 
(1 — Re detR) = 0 as a new Q Ward identity when G > 0. Across a wide range of Ai at on 4 4 
lattices, the improved action leads to smaller (1 — Re detR) than the unimproved action produces. 
Lines connect points with fixed n or G (cf. table 1). While larger values of n or G trivially 
suppress (1 — Re detR), we want to use the smallest acceptable values to stay as close as possible 
to undeformed Af = 4 SYM. 


range of Ai at on small 4 4 lattices. While the improved action produces smaller values for 
this quantity, on this lattice volume the results from the two actions are not dramatically 
different, especially at weaker couplings. From figure 4 and table 1 we can also see that 
larger values of the auxiliary coupling G trivially reduce (1 — Re det V). However, we still 
want to minimize G to remain as close as possible to undeformed Af = 4 SYM. 

The real contrast between 5 un i mp and Si mp appears when we consider the ( a/L ) —> 0 
continuum limit in figure 5. As for the bosonic action deviations, the improved action 
produces (1 — Re det V) oc (a/L) 2 for all investigated’t Hooft couplings 1 < Ai at < 4. The 
specific exponents p resulting from power-law fits oc ( a/L) p are collected in table 2. For 
the unimproved action with Ai at = 1 it is not even clear that the U(l) sector is actually 
decoupling as L increases within the range 4 < L < 12. The corresponding power-law fit 
produces an exponent p = 0.0027(21) that is barely distinguishable from zero. 

The last Ward identity we consider involves the ijip a fermion bilinear and was intro¬ 
duced in ref. [9]: ( QO) = 0 where O = Tr and so 


QO = Tr 


dYUaU a 


-Tr 


7?y ^ Ip gif c 


= D — F. 


(4.1) 


We use the shorthand “D” and “F" because the first term depends on the modified equa¬ 
tions of motion for the auxiliary field d, eq. 3.3, while the second involves the fermion bilin¬ 
ear. In figure 6 we show the usual power-law continuum extrapolations of the violations of 
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Figure 5. Power-law continuum extrapolations of (1 — Re defP) on logarithmic axes for the 
unimproved action with Ai a t = 1 and the improved action with 1 < Ai a t < 4. For Sj mp this quantity 
vanishes ~(a/L) 2 for each coupling (table 2), as expected when the U(l) sector decouples from the 
SU(iV) target theory. The gap between Ai at = 1 and the stronger couplings is due to the larger 
values of p that had to be used for Ai a t > 2 (cf. table 1). For 5 U nimp the corresponding exponent 
p = 0.0027(21) is not significantly non-zero, implying that L > 12 is needed to reliably observe 
unimproved simulations approaching the SU(fV) continuum theory. 



Figure 6. Power-law continuum extrapolations of normalized ry?/< a bilinear Ward identity violations 
(eq. 4.1) on logarithmic axes for the unimproved action with Ai a t = 1 and the improved action with 
1 < Ai a t < 4. This quantity vanishes oc ( a/L) p in the ( a/L ) —>• 0 continuum limit, withp = 0.040(17) 
for the unimproved action. For the improved action larger violations on small lattice volumes are 
quickly compensated by the much larger p«2 collected in table 2. The gap between Ai a t = 1 and 
the stronger couplings is due to the larger values of fi that had to be used for Ai at > 2 (cf. table 1). 
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| (s B ) ~ 9N 2 /2\ 

(1 — Re det V) 2 

/ QO 

9N 2 /2 

\D + F 

Unimproved, Ai at = 1 

0.015(4) 

0.0027(21) 

0.040(17) 

Improved, Ai at = 1 

2.42(13) 

1.94(15) 

2.21(15) 

Improved, Ai at = 2 

1.99(9) 

1.93(6) 

2.20(7) 

Improved, Ai at = 3 

1.99(9) 

2 .10(6) 

2.46(7) 

Improved, A lat = 4 

1.90(6) 

1.97(6) 

2.32(6) 


Table 2. Results for the exponents p obtained by fitting the listed lattice artifacts to power-law 
continuum extrapolations oc ( a/L) p as shown in figures 2, 3, 5 and 6. All exponents for the improved 
action are p ~ 2, suggesting that we have achieved 0{a ) improvement as discussed in the text and in 
appendix B. Exponents for the unimproved action are much smaller, especially for (1 — Re detV ), 
which is only protected by a Q Ward identity for Si mp with non-zero G. 
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Figure 7. Comparison of the D and F terms in eq. 4.1 for the T]i/j a bilinear Ward identity at fixed 
Aiat = 1. For the unimproved action both terms are fairly insensitive to the volume, with comparable 
uncertainties. The fermion bilinear F term behaves similarly for the improved action, while the pure- 
gauge D term changes much more significantly and ends up with much larger jackknife uncertainties 
due to the modified equations of motion for the auxiliary field d (eq. 3.3). 


this Ward identity, normalized by + F). As for the other quantities considered above, 
the improved action produces Ward identity violations that vanish ~(a/L) 2 in the contin¬ 
uum limit, while the exponent for the unimproved action is much smaller, p = 0.040(17). 
Precise values of p for the improved action at each 1 < Ai a t < 4 are collected in table 2. 

One curious feature of figure 6 is that the improved action leads to much larger ( QO) 
on small lattice volumes, especially 4 4 . Although this is quickly compensated by the more 
rapid approach to the continuum limit as L increases, it is a striking change compared to 
figures 2 and 5. Contrary to our expectations, we find that the fermion bilinear term F 
remains fairly insensitive to the volume for the improved action. It is the pure-gauge term D 
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that now changes significantly as L increases, steadily approaching F as shown in figure 7. 
The improved action also leads to large fluctuations in the D term, producing jackknife 
uncertainties more than an order of magnitude larger than those of F (even though we 
use only three stochastic sources to measure F on each saved gauge configuration). Not 
surprisingly given these results, the fermion BCs have no statistically significant effect on 
these Ward identity violations, which change from 0.056(2) to 0.063(9) when we switch 
from anti-periodic to periodic BCs for (Ai a t,/x, G) = (1,0.4,0.1) on 4 4 lattices. 



| (s B ) - 9iV 2 /2| 

(1 — Re det V) 

2 

/ 2° \ 



9N 2 /2 

\D + F/ 


U(2) 

0.0062(11) 

0.0212(49) 

0.0849(176) 

U(3) 

0.0028(5) 

0.0060(19) 

0.0254(57) 

U(4) 

0.0024(2) 

0.0036(10) 

0.0162(31) 


Table 3. Bosonic action deviations, plaquette determinant fluctuations and r/ipa bilinear Ward 
identity violations from the improved action for gauge groups U(1V) with N = 2, 3 and 4 on 4 4 
lattices with fixed (Ai at ,/r,G) = (1,0.4,0.05). The results are reasonably consistent with 1/N 2 sup¬ 
pression; some deviations from such scaling are attributable to the anti-periodic fermion temporal 
BCs. 


Finally, in ref. [9] we reported that Ward identity violations for the unimproved action 
decrease oc 1/N 2 for gauge group U(1V). While all of the results discussed above are for 
U(2), our initial explorations of larger N indicate that the improved action also exhibits 
this 1 /N 2 suppression of supersymmetry breaking. In table 3 we collect the bosonic action 
deviations, (1 — Re det V ). and r)'i/> a bilinear Ward identity violations from improved 4 4 
calculations with N = 2 , 3 and 4 at fixed (Ai a t, yu, G) = (1,0.4,0.05). These results are 
reasonably consistent with 1/N' 2 scaling, though some of them become so small that they 
appear sensitive to the anti-periodic fermion temporal BCs. 

5 Conclusion 

In this paper we have presented a procedure for lifting the U(l) flat directions in lattice 
J\T = 4 SYM in a way that preserves the exact lattice supersymmetry. These U(l) flat 
directions are a necessary feature of the supersymmetric lattice construction, and must be 
regulated to ensure that the path integral is well defined and for the lattice theory to have 
the correct naive continuum limit. Our procedure involves a deformation of the moduli 
space of the lattice theory, based on modifying the equations of motion of an auxiliary 
field. 

We applied this procedure to construct a new lattice action that we compared against 
the “unimproved” action used in earlier work. We saw that the new action reduces viola¬ 
tions of supersymmetric Ward identities by up to an order of magnitude on small L 4 lattice 
volumes, and more importantly we found that these violations now fall rapidly towards zero 
oc 1/L 2 as L increases. Indeed, we argue that the new action is effectively 0(a)-improved 
since there are no dimension-5 operators compatible with the lattice symmetries. These 
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numerical benefits appear to be well worth the roughly doubled computational expense, 
and we recommend that the improved action introduced in this work be used in future 
studies of lattice Af = 4 SYM. 
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Appendix A: Another possible lattice action and its difficulties 

As discussed at the end of section 3.2, it is enticing to consider the possibility of con¬ 
structing a lattice action that includes both of the necessary scalar potential and plaquette 
determinant deformations while still maintaining exact Q supersymmetry. This requires 
making each term non-negative so that it cannot cancel out the other, leading us to consider 

r]{n) (v[ ^U a (n) + G^2\det V a b(n) - l\ 2 l N (A.l) 

V a^b 

+ B 2 J2 [U a (n)U a (n)\ - l) I N 

n. ' ' 


with two tunable parameters B and G. Although we use the same symbol G for its coupling, 
the plaquette determinant term has changed compared to eq. 3.6. 

The modified equations of motion for the auxiliary field d are now 



(A.2) 


which produce the Q Ward identity 

dEE |det V a b(n) - 1| 2 \ + B 2 (ee O[«.(»*.(»>] -Oh (a.3) 

\ n a^b / \ n a ' / / 

from eq. 3.4. Every term in these sums is non-negative, so the only way for this Ward 
identity to be satisfied is for each term to vanish independently on each lattice site. That 
is, even though the lattice action is now fully Q-invariant (with periodic BCs), reaching 
a supersymmetric vacuum requires simultaneously satisfying at least 15V non-trivial con¬ 
straints, where V is the lattice volume. We will soon see that this is not feasible in practice, 



Q Tr 


77(77.) 


V[ ] U n 


(n 


Q Tr 
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which leads us to name this possibility ‘the over-constrained action’. Written in full, it is 


Saver — S exac (. + S c l ose( j 

Q/f _ Q/ Q 

^ exact ^det *^pot 

+—E Tr 


(A.4) 


2A l 


at 


F a b{n)F ab {n) - Xab(n)V^ip b] {n) - r]{n)V\ > a (n) 




+ \[ v< a ] Ua(n) + GJ2 Idet P ab (n) - l| 2 Ihv + B 2 ^ (j-Tr [U a (n)U a (n)\ - l] I N 

\ a^b a, ' ' . 


ql 

“-Met 


V 7 —2G Tr [rj(n)j Y] det P a 6 (n) [det P ba (n) - 1] 

^Vat ^ 

X Tr [U^ 1 (n^bin) + U^ l (n+ '{l b )^ a {n+ ^ b )] 


Spot — 


N 2 B 2 


2Ai at N 

with S c i OS ed still given by eq. 2.3. 


E Tr W 74 )] E ( [Ua(n)U a (n)] - 1 ) Tr [^ a (n)U a {n)\ 


Over-constrained action 


Volume 

Alat 

B 

G 

4 4 

0.5 

0.5 

0.10 

4 4 

1.0 

0.4 

0.05 

4 4 

2.0 

0.5 

0.05 

4 4 

3.0 

0.5 

0.10 

4 4 

4.0 

0.5 

0.15 

4 4 

5.0 

0.5 

0.20 

4 4 

6.0 

0.5 

0.25 

4 4 

7.0 

0.5 

0.25 

4 4 

8.0 

0.5 

0.30 

4 4 

9.0 

0.5 

0.30 

4 4 

10.0 

0.5 

0.35 

6 4 

1.0 

0.3 

0.05 

8 4 

1.0 

0.2 

0.05 

12 4 

1.0 

0.15 

0.05 


Table 4. Auxiliary couplings used to obtain the results for the over-constrained action in figures 8- 


With the over-constrained action defined, we can carry out the same parameter space 
exploration and investigation of the ( a/L ) —> 0 continuum limit as for the improved action 
in section 4. The resulting lattice ensembles are listed in table 4 and produce the results 
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^lat 

Figure 8. Deviations of the bosonic action from its exact supersymmetric value for the over¬ 
constrained action in addition to the other results in figure 1 for a wide range of’t Hooft coupling 
Ai a t on 4 4 lattices. Although S ovel leads to smaller deviations than the unimproved action produces, 
it does not behave as well as the improved action. Lines connect points with fixed B (cf. table 4). 


shown in figures 8-13. Since these figures present the same quantities that we discussed at 
some length in section 4, here we focus only on the new results from the over-constrained 
action. First, the scan of bosonic action deviations across a wide range of’t Hooft coupling 
Ai a t on small 4 4 lattices in figure 8 reveals that S over still produces non-zero Q Ward identity 
violations despite avoiding all explicit supersymmetry breaking in the action. In fact, S over 
behaves worse than the improved action, though it still provides some benefit compared 
to 5 i in im p. Of course, the anti-periodic temporal BCs for the fermions do introduce some 
explicit Q breaking, but as for the other actions this is responsible for only a small portion 
of the deviations shown. For B,G) = (1,0.5, 0.1) on 3 3 x 4 lattices, switching from 
anti-periodic to periodic BCs only reduces the bosonic action deviation from 0.0260(11) to 
0.0224(5). 

The picture becomes worse in figure 9, where we consider the ( a/L ) —> 0 continuum 
extrapolation of the bosonic action deviations produced by the over-constrained action, 
and find that they vanish oc ( a/L) v with p = 0.14(4). While this exponent is an order of 
magnitude larger than that for the unimproved action (table 2 ), it is a far cry from the p ~ 2 
obtained for the improved action. That is, despite preserving the exact Q supersymmetry 
at the level of the action, 5 over is not effectively O(a) improved. 

The situation is similar when we consider (1 — Re det V), again averaging the plaquette 
determinant over all orientations and lattice sites even though the Ward identity in eq. A.3 
should protect det V a b(n) ~ 1 site by site. In figure 10 we see that the over-constrained 
action often does a worse job controlling deCP than the unimproved action, even though 
we end up having to use relatively large G as listed in table 4. We are also unable to 
observe the U(l) sector decoupling as we approach the continuum limit in figure 11, where 
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Figure 9. Power-law continuum extrapolations of bosonic action deviations on logarithmic axes for 
the over-constrained action in addition to the other results in figure 2 at fixed Aj a t = 1. Although 
S'over produces L = 4 deviations comparable to those of the improved action, as L increases these 
only decrease oc ( a/L) p with a much smaller power p = 0.14(4). 



^lat 


Figure 10. (1 — Re defP) for the over-constrained action in addition to the other results in 

figure 4 for a wide range of’t Hooft coupling Ai at on 4 4 lattices. Even though this quantity should 
be protected by a Q Ward identity (eq. A. 3), S ovel often behaves worse than the unimproved action 
despite using relatively large G (cf. table 4). Lines connect points with fixed G. 


fitting (1 — Re defP) oc ( a/L) p produces p = 0.0021(43) indistinguishable from zero. Since 
(1 — Re detR) should be protected by a Q Ward identity, this may also raise concerns 
about the restoration of the Q supersymmetry from S ovei in the continuum limit. 

Finally, we consider the r/?/> a bilinear Ward identity violations in figures 12 and 13. The 
over-constrained action clearly produces the worst results for this quantity, apparently due 
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Figure 11. Power-law continuum extrapolations of (1 — Re detP) on logarithmic axes for the 
over-constrained action in addition to the other results in figure 5. Not only does S ovel produce 
larger plaquette determinant fluctuations even though this quantity should be protected by a Q 
Ward identity (eq. A.3), the exponent p = 0.0021(43) from fitting (1 — Re detP) oc ( a/L) p is 
consistent with zero. 



Figure 12. Power-law continuum extrapolations of normalized r)il) a bilinear Ward identity viola¬ 
tions (eq. 4.1) on logarithmic axes for the over-constrained action in addition to the other results 
in figure 6. The modified equations of motion for d in eq. A.2 produce very large QO = D — F that 
vanish very slowly in the continuum limit, oc ( a/L) p with p = 0.031(1). 


to the significant modifications of the auxiliary held equations of motion in eq. A.2. These 
modifications lead to extremely large values for the pure-gauge term D in eq. 4.1, while 
the fermion bilinear term F from S ovei is not changed so significantly. Although D does 
approach F as the volume increases, the power-law fit oc ( a/L) p in figure 12 still produces 
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Figure 13. Comparison of the D and F terms in the r/i/a bilinear Ward identity (eq. 4.1) for the 
over-constrained action in addition to the other results in figure 7 at fixed Ai a t = 1. The modified 
equations of motion for d in eq. A.2 produce pure-gauge terms D that are now much larger than 
the fermion bilinear terms F , and only slowly decrease as L increases. 


a very small p = 0.031(1), worse even than the unimproved action. 

In summary, despite maintaining exact Q supersymmetry the over-constrained action 
produces much more severe Ward identity violations than the improved action, and results 
in a much slower approach to the {a/L) —> 0 continuum limit. The reason for this behavior 
seems to be that S over requires the system to satisfy the 1517 non-trivial constraints in 
eq. A.3 in order to reach a supersymmetric vacuum, which is not possible in practice. 
Since the Ward identities are violated even with fully periodic BCs, 5 ove r appears to result 
in spontaneous supersymmetry breaking, which is forbidden in continuum M = 4 SYM 
where the Witten index is non-zero. This suggests that the over-constrained lattice system 
is changing the Witten index by introducing new lattice states that don’t exist in the 
continuum, which may be related to the slow approach to the continuum limit seen in 
figures 9, 11 and 12. 

While it may be worthwhile to continue searching for ways to further reduce or elim¬ 
inate soft supersymmetry breaking in lattice J\f = 4 SYM, the alternate lattice action 
considered in this appendix is not viable. In the context of practical numerical computa¬ 
tions, the improved action introduced in eq. 3.10 appears hard to beat. 

Appendix B: Dimension-5 operators in lattice J\f = 4 SYM 

In this appendix we extend the enumeration of renormalizable operators in lattice A/" = 4 
SYM from ref. [10] to consider dimension-5 operators as well . 3 In addition to respecting 
the Q supersymmetry, lattice gauge invariance and the £5 point group symmetry of the A\ 
lattice, these operators must also be uncharged under an additional U(l) “ghost number” 

J We thank Joel Giedt for discussions on this subject. 
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symmetry [1, 10]. This ghost number symmetry is hidden in the notation used to this 
point, and only becomes visible when the lattice fields are decomposed into irreducible 
representations of £5 by means of the 5x5 orthogonal matrix P discussed in refs. [3, 9]. 
Schematically, 

U a 4 (B.l) 

Xab ^ (Xnv^n) 

where fi = 0 • ■ ■ 3. These fields have the ghost charges 

4> 0 r) Xvv Vv V 

0 2 0 -2 1 -1 1 -1 1 . 

Now we can begin by confirming that these symmetries permit no dinrension-5 Q-closed 
operators. Generic Q-closed terms involve only the fields that cannot be obtained through 
a supersymmetry transformation, namely 77, U a and Xab■ Two dimension-5 combinations 
of these fields are invariant under the £5 point group symmetry: 

XabV^V^rj €abode (B.2) 

The first of these terms vanishes because Xab is anti-symmetric under the interchange a O b 
while T>^T>^ is symmetric. The e tensor in the other term requires that (j> appear, leading 
to a non-zero ghost charge that forbids the operator. 

Dinrension-5 Q-exact operators can take the form Q TV [\E , i\E , 2 ’I , 3 ] in addition to the 
Q TV f(U a ,U a ,d)] and Q {TV [ 77 ] TV [f(U a ,U a ,d)]} considered in ref. [10]. Here each T 
stands for one of the fermion fields 77 , -0a or Xab and f(U a ,U a ,d ) has to be a dimension- 
3 combination of the bosonic fields. It is easy to see that there are no dinrension-5 
operators of the form Q {TV [ 77 ] TV [f(U a ,U a , d)] }. Such terms require an S' 5 -invariant 
f(U a ,U a ,d) —> f(U a U a ,d), which can only have dim[/] = 2 , 4 or higher. While this also 
rules out Q TV [T !f(JA a ,U a ,d)\ when T = 77 , the £5 symmetry allows terms with T = ij) a 
and T = Xab- 

Q TV [t/> a U a d] Q Tr [^ a U a U b U b ] Q Tr [e abcde XalUMdUe] ■ (B.3) 

The first term vanishes since Q annihilates each of 7 /> a , U a and d (eq. 2.5), while the e tensor 
in the last term requires that 4> appear, leading to a non-zero ghost charge. The middle 
term also has a non-zero ghost charge, and was already discussed in ref. [10]. Finally, the 
only £ 5 -invariant Q Tr [T'iT 2 v I / 3 ] terms are 

Q Tr [777777] Q Tr [i’aAXab] , (B.4) 

both of which are forbidden by ghost charge conservation. 

So we conclude that exact Q supersymmetry and the other symmetries of lattice J\T = 4 
SYM forbid all dimension-5 operators. When supersymmetry breaking is sufficiently small, 
we may therefore expect the system to be effectively 0{a) improved, explaining the behavior 
of the improved action observed in section 4. 
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Appendix C: New results for the pfaffian phase 


The improved and over-constrained lattice actions both add the plaquette determinant to 
the fermion operator V in such a way as to break the —> r) + cl^v shift symmetry and lift 
the corresponding U(l) fermion zero mode. Cf. S^et i n eq. 3.10 and S' det in eq. A.4. This 
has the serendipitous effect of allowing us to compute the pfaffian of V with fully periodic 
BCs, which is not possible for the formal or unimproved lattice actions. 

As discussed in more detail by ref. [9], pfV = |pfP|e* a results from integrating over 
the fermion fields in the lattice path integral and is potentially complex, a/0. We carry 
out phase-quenched computations that determine observables (O) without including the 
e*“ factor in the path integral, which we can then reintroduce through phase reweighting, 


(O) 


(Oe ia ) 

\ 'pq 


(C.l) 


In ref. [9] we found that e*“ ~ 1 with anti-periodic BCs, which makes the phase reweighting 
step irrelevant. If the phase had fluctuated to the extent that (e ia ) pq became consistent 
with zero, then we would have faced a sign problem. 

Because we used the unimproved action in ref. [9] , we were not able to investigate the 
pfaffian phase with periodic BCs. In figure 14 we compare the pfaffian phase with periodic 
vs. anti-periodic BCs for both the improved action and the over-constrained action. In both 
cases anti-periodic BCs produce e ia ~ 1 as we found for the unimproved action. Taking 


N — 4 SYM, U(2) IAnti-periodic BCs + 
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Figure 14. Scatter plots of pfaffian phase measurements from independent 3 3 x4 phase-quenched en¬ 
sembles using either periodic or anti-periodic temporal BCs for the fermion fields, for the improved 
action with (Ai at ,/z,G) = (1,0.8,0.1) (left) and the over-constrained action with (Ai at ,.B,G) = 
(1,0.5, 0.1) (right). For both actions anti-periodic BCs produce e za ~ 1 as we found for the unim¬ 
proved action in ref. [9]. Periodic BCs, however, lead to uncontrolled fluctuations and ( e *“) p9 
consistent with zero, indicating a sign problem. 
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the ensemble average, 


(e ia ) pq = 0.9977(4) + *0.0003(82) for S imp 
(e ia ) pq = 0.9952(9) + *0.0020(118) for 5 over . 

With periodic BCs, however, the pfaffian phase fluctuates enough to average to zero for 
both actions, indicating a sign problem: 

(■e ia ) pq = 0.071(82) - *0.029(88) for 5 imp 
(e ia ) pq = -0.041(80) + *0.086(78) for 5 over . 

It is not yet clear to us why the pfaffian is so sensitive to the temporal fermion BCs. 
Even more mysteriously, all other observables change very little between the ensembles 
generated with periodic or anti-periodic BCs, despite the apparent presence of a sign prob¬ 
lem in the former case and its absence in the latter. For example, the bosonic action 
deviations from the 5j mp ensembles that produce the left plot in figure 14 are 0.0169(8) 
with anti-periodic BCs and 0.0115(8) with periodic BCs, just what we would expect in the 
absence of a sign problem. This numerical evidence suggests that lattice A/" = 4 SYM does 
not suffer from a sign problem even when {e ia ) pq is consistent with zero. We continue to 
search for ways to analytically understand this unexpected result. 
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